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We combine results from CDF and DO's direct searches for the standard model (SM) Higgs boson 
(H) produced in pp collisions at the Fermilab Tevatron at ^/s — 1.96 TeV, focusing on the decay 
H — >■ 77. We compute upper limits on the Higgs boson production cross section times the decay 
branching fraction in the range 100 < mn < 150 GeV/c^, and we interpret the results in the 
context of the standard model. We use the MSTW08 parton distribution functions and the latest 
theoretical cross section predictions when testing for the presence of a SM Higgs boson. With 
^L.' . datasets corresponding to 7.0 fb~^ (CDF) and 8.2 fb~^ (DO), the 95% C.L. upper limits on Higgs 

• ' boson production is a factor of 10.5 times the SM cross section for a Higgs boson mass of 115 GeV/c^. 

o 



Preliminary Results 



The Tevatron New-Phenomena and Higgs Working Group can be contacted at TEVNPHWG@fnal.gov. More information can be found 
at http://tevnphwg.fnal.gov/ 



I. INTRODUCTION 

The search for a mechanisni for electroweak symmetry breaking, and in particular for a standard model (SM) Higgs 
boson, has been a major goal of particle physics for many years, and is a central part of the Fermilab Tevatron physics 
program, and also that of the Large Hadron Collider. Recently, the ATLAS and CMS Collaborations at the Large 
Hadron Collider have released results on searches for the standard model Higgs boson decaying to W^W~ [ll, 0], 
ZZ ^] , and 77 [3, |5| . Both the CDF and DO Collaborations have performed new combinations [6|-|8[ of multiple direct 
searches for the SM Higgs boson. The sensitivities of these new combinations significantly exceed those of previous 
combinations [9lllQ'|. 

In this note, we combine the most recent results of searches for the SM Higgs boson produced in pp collisions 
at ^/s = 1.96 TeV, where the decay of the Higgs boson is restricted to pairs of photons: H -^ 77. We consider all 
four SM production modes, gluon-gluon fusion: gg — >■ H, associated production of a Higgs boson with a W boson: 
qq — >■ WH, associated production of a Higgs boson with a Z boson: qq — » ZH, and vector-boson fusion: qq -^ q'q'H. 
The CDF search 11 [ is performed with 7.0 fb~^ and the DO search [l^] is performed with 8.2 fb^^ of collision data. 



The data for the CDF search are divided into four categories: Central-Central (CC), Central-Plug (CP), Central- 
Central Conversion (CC Conv), and Central-Plug conversion (CP Conv), which are based on whether the photon(s) 
are detected in the central electromagnetic (EM) calorimeter {\rj\ < 1.1) or the plug EM calorimeter (1.2 < \ri\ < 2.8), 
and whether or not one of the central photons is detected via its conversion products e+e". The reconstructed 
mass resolution is approximately 2.8 GeV/c^ for the CDF channels. The acceptances for Higgs boson decay are 
approximately 13% for the Central-Central category, 16% for the Central-Plug category, 2.9% for the Central-Central 
conversion category, and 1.8% for the Central-Plug Conversion category, varying by ±10% (relative) depending on the 
mass of the Higgs boson and the production mechanism. We assume the SM mixture of the production mechanisms. 
The discriminant variable is the reconstructed candidate mass m^-y, and the background is parameterized as a smooth 
function of m^^^ which is fit to the data outside of a window centered on the signal under test. 

The DO search requires events with at least two photon candidates with Et > 25 GeV and |?7| < 1.1. The diphoton 
mass resolution is ~ 3 GeV/c^. The contribution of jets misidentified as photons is reduced by combining information 
sensitive to differences in the energy deposition from these particles in the tracker, calorimeter and central preshower 
in a neural network. In this latest iteration of the analysis boosted decision trees, rather than the diphoton invariant 
mass, are used as the final discriminating variable. The transverse energies of the leading two photons along with 
azimuthal opening angle between them and the di-photon invariant mass and transverse momentum arc used as input 
variables. A sizeable improvement in sensitivity beyond that achieved with the invariant mass is obtained. 

II. SIGNAL PREDICTIONS AND UNCERTAINTIES 

We normalize our Higgs boson signal predictions to the most recent highest-order calculations available, for all 
production processes considered. The largest production cross section, a{gg — >■ H), is calculated at next-to-next-to- 
leading order (NNLO) in QCD with soft gluon resummation to next-to-next-to-leading-log (NNLL) accuracy, and 
also includes two- loop electroweak effects and handling of the running b quark mass |13l . Il4| . The numerical values 
in Table U are updates |15| of these predictions with rrit set to 173.1 GeV/c^ ^i?], and an exact treatment of the 
massive top and bottom loop corrections up to next-to-leading-order (NLO) + next-to-leading-log (NLL) accuracy. 
The factorization and renormalization scale choice for this calculation is fj^p = U r = J't-h - These calculations are 
refinements of earlier NNLO calculations of the gg ^ H production cross section J17H19| . Electroweak corrections 
were computed in Refs. |20l.l2l|. Soft gluon resummation was introduced in the prediction of the gg ^ H production 
cross section in Ref. |23 |. 

The gg ^ H production cross section depends strongly on the gluon parton density function, and the accompanying 
value oiasiq"^). The cross sections used here are calculated with the MSTW08 NNLO PDF set [2^, as recommended by 
the PDF4LHC working group \2i\. We follow the PDF4LHC working group's prescription to evaluate the uncertainties 
on the gg —i' H production cross section due to the PDFs. This prescription is to evaluate the predictions of (j{gg — >■ H) 
at NLO using the global NLO PDF sets CTEQ6.6 [H, MSTW08 [13, and NNPDF2.0 [1^, and to take the envelope 
of the predictions and their uncertainties due to PDE-l-a^ for the three sets as the uncertainty range at NLO. The 



TABLE I: The production cross sections and decay branching fractions for the SM Higgs boson assumed for the combination. 
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100 


1821.8 


291.90 


169.8 


100.1 


0.159 


105 


1584.7 


248.40 
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0.178 


110 


1385.0 


212.00 


125.7 


85.1 


0.197 


115 


1215.9 


174.50 


103.9 


78.6 


0.213 


120 


1072.3 


150.10 


90.2 


72.7 


0.225 


125 


949.3 


129.50 


78.5 


67.1 


0.230 


130 


842.9 


112.00 


68.5 


62.1 


0.226 
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750.8 


97.20 


60.0 


57.5 


0.214 


140 


670.6 


84.60 


52.7 


53.2 


0.194 


145 


600.6 


73.70 


46.3 


49.4 


0.168 


150 


539.1 


64.40 


40.8 


45.8 


0.137 



ratio of the NLO uncertainty range to that of just MSTW08 is then used to scale the NNLO MSTW08 PDF+as 
uncertainty, to estimate a larger uncertainty at NNLO. This procedure roughly doubles the PDF+as uncertainty from 
MSTW08 at NNLO alone. 

We also include uncertainties on cr^gg — >■ H) due to uncalculated higher order processes by following the standard 
procedure of varying the factorization renormalization scales up and down together by a factor k = 2. We treat the 
scale uncertainties as 100% correlated between jet categories and between CDF and DO, and also treat the EDF+a^ 
uncertainties in the cross section as correlated between jet categories and between CDF and DO. We treat however 
the PDF+tts uncertainty as uncorrelated with the scale uncertainty. 

We include all significant Higgs boson production modes in our searches. Besides gluon-gluon fusion through virtual 
quark loops, we include Higgs boson production in association with a W or Z vector boson, and vector boson fusion 
(VBF). We use the WH and ZH production cross sections computed at NNLO in Ref. [27| . This calculation starts 
with the NLO calculation of v2hv [28| and includes NNLO QCD contributions [29], as well as one- loop electroweak 
corrections [30] ■ We use the VBF cross section computed at NNLO in QCD in Ref. [31|. Electroweak corrections to 
the VBF production cross section are computed with the HAWK program [321 ' ^^^ are small (3% or less) for the Higgs 
boson mass range considered here. The VBF cross sections in Table U do not include the electroweak corrections. 

In order to predict the distributions of the kinematics of Higgs boson signal events, CDF and DO use the pythia [Sj] 
Monte Carlo program, with CTEQ5L and CTEQ6L [34| leading-order (LO) parton distribution functions. 

The Higgs boson decay branching fraction Br{H — )■ 77) used here is obtained from the Handbook [35]. The 
branching fractions are computed by evaluating the partial decay widths of on-shell Higgs bosons to all possible final 
states allowed in the SM, and evaluating the fractions of the total decay width for each process. For all processes 
except H — >■ W^W~ and H — >• ZZ, HDECAY ^ is used to compute the partial widths, hdecay includes relevant 
higher-order QCD corrections to decays into quarks and gluons, and NLO electroweak corrections are included in the 
i/ — > 77 and H ^ gg processes. For H — > W^W~ and H -^ ZZ, the Monte Carlo generator Prophecy4f [37| is 
used. It computes partial widths for Higgs boson decays at NLO, including NLO QCD and electroweak corrections 
and all interferences, for four-fermion final states of Higgs boson decay. More details are available in Ref. [33 • The 
branching fraction Br{H —> 77) is listed in Table H] as a function of ruH- 

III. COMBINING CHANNELS 

To gain confidence that the final result does not depend on the details of the statistical formulation, we perform 
two types of combinations, using Bayesian and Modified Frequentist approaches, which yield limits on the Higgs 
boson production rate that agree within 5% at each value of mn, and within 1% on average. Both methods rely on 
distributions in the final discriminants, and not just on their single integrated values. Systematic uncertainties enter 



on the predicted number of signal and background events as well as on the distribution of the discriminants in each 
analysis ("shape uncertainties"). Both methods use likelihood calculations based on Poisson probabilities. 

Both methods treat the systematic uncertainties in a Bayesian fashion, assigning a prior distribution to each 
source of uncertainty, parameterized by a nuisance parameter, and propagating the impacts of varying each nuisance 
parameter to the signal and background predictions, with all correlations included. A single nuisance parameter 
may affect the signal and background predictions in many bins of many different analyses. Independent nuisance 
parameters are allowed to vary separately within their prior distributions. Both methods use the data to constrain 
the values of the nuisance parameters, one by integration, the other by fitting. These methods reduce the impact 
of prior uncertainty in the nuisance parameters thus improving the sensitivity. Because of these constraints to the 
data, it is important to evaluate the uncertainties and correlations properly, and to allow independent parameters to 
vary separately, otherwise a fit may overconstrain a parameter and extrapolate its use improperly. The impacts of 
correlated uncertainties add together linearly on a particular prediction, while those of uncorrelated uncertainties are 
convoluted together, which is similar to adding in quadrature. 

A. Bayesian Method 

Because there is no experimental information on the production cross section for the Higgs boson, in the Bayesian 
technique Q [331 ^^ assign a flat prior for the total number of selected Higgs events. For a given Higgs boson mass, the 
combined likelihood is a product of likelihoods for the individual channels, each of which is a product over histogram 
bins: 

Nc Nt «np 

c{R, s, b\n, e) X ^{9) = n n ^^l' ^"''" /"^^- • "" n ^"''^' (i) 

i=lj=l k=l 

where the first product is over the number of channels (A^c), and the second product is over N^ histogram bins 
containing rnj events, binned in ranges of the final discriminants used for individual analyses, such as the reconstructed 
dijet mass or neural-network outputs. The parameters that contribute to the expected bin contents are ^^ = 
R X Sij{9) + hij{9) for the channel i and the histogram bin j, where Sy and bij represent the expected signal and 
background in the bin, and i? is a scaling factor applied to the signal to test the sensitivity level of the experiment. 
Truncated Gaussian priors are used for each of the n^p nuisance parameters 9k, which define the sensitivity of the 
predicted signal and background estimates to systematic uncertainties. These can take the form of uncertainties on 
overall rates, as well as the shapes of the distributions used for combination. These systematic uncertainties can be 
far larger than the expected SM Higgs boson signal, and are therefore important in the calculation of limits. The 
truncation is applied so that no prediction of any signal or background in any bin is negative. The posterior density 
function is then integrated over all parameters (including correlations) except for R, and a 95% credibility level upper 
limit on R is estimated by calculating the value of R that corresponds to 95% of the area of the resulting distribution. 

B. Modified Frequentist Method 

The Modified Frequentist technique relies on the CLg method, using a log-likelihood ratio (LLR) as test statistic [21 : 

ii,?^_21n4^^^, (2) 

p(data|iJo) 

where Hi denotes the test hypothesis, which admits the presence of SM backgrounds and a Higgs boson signal, while 
Hq is the null hypothesis, for only SM backgrounds. The probabilities p are computed using the best-fit values of the 
nuisance parameters for each pseudo-experiment, separately for each of the two hypotheses, and include the Poisson 
probabilities of observing the data multiplied by Gaussian priors for the values of the nuisance parameters. This 



technique extends the LEP procedure [38| which does not involve a fit, in order to yield better sensitivity when 
expected signals are small and systematic uncertainties on backgrounds are large [39j. 

The CLg technique involves computing two p- values, CLg+b and CLb- The latter is defined by 

1 - CLb - p{LLR < iLi?obs|i?o), (3) 

where LLRohs is the value of the test statistic computed for the data. 1 — CLb is the probability of observing a 
signal-plus-background-like outcome without the presence of signal, i.e. the probability that an upward fluctuation of 
the background provides a signal-plus-background-like response as observed in data. The other p-value is defined by 

CLs+b = piLLR > LLR,i,,\Hi), (4) 

and this corresponds to the probability of a downward fluctuation of the sum of signal and background in the data. 
A small value of CLg+b reflects inconsistency with Hi. It is also possible to have a downward fluctuation in data 
even in the absence of any signal, and a small value of CLg+b is possible even if the expected signal is so small 
that it cannot be tested with the experiment. To minimize the possibility of excluding a signal to which there is 
insufficient sensitivity (an outcome expected 5% of the time at the 95% C.L., for full coverage), we use the quantity 
CLs = CLs+b/CLb. If CLs < 0.05 for a particular choice of Hi, that hypothesis is deemed to be excluded at the 
95% C.L. In an analogous way, the expected CLb, CLs+b and CLg values are computed from the median of the LLR 
distribution for the background-only hypothesis. 

Systematic uncertainties are included by fluctuating the predictions for signal and background rates in each bin of 
each histogram in a correlated way when generating the pseudo-experiments used to compute CLg+b and CLb. 

C. Systematic Uncertainties 

Systematic uncertainties differ between experiments and analyses, and affect the rates and shapes of the predicted 
signal and background in correlated ways. The combined results incorporate the sensitivity of predictions to values 
of nuisance parameters, and include correlations between rates and shapes, between signals and backgrounds, and 
between channels within experiments and between experiments. More discussion on this topic can be found in the 
individual analysis notes |11| and [l2|. Here we describe only the largest contributions and correlations between and 
within the two experiments. 

1. Correlated Systematics between CDF and DO 

The uncertainties on the measurements of the integrated luminosities are 6% (CDF) and 6.1% (DO). Of these 
values, 4% arises from the uncertainty on the inelastic pp scattering cross section, which is correlated between CDF 
and DO. 

In both CDF and DO's iJ ^- 77 analyses, the dominant background yields are calibrated with data control samples, 
and thus the systematic uncertainties on the background rates and shapes are considered uncorrelated between the 
two collaborations' results. 

The theoretical cross section uncertainties on the inclusive gg -^ H production cross section (^ 14%), the WH 
and ZH cross sections (~ 6%), and the VBF cross section ('^ 5%) are shared between CDF and DO and are used 
in the computation of the limits scaled to the SM prediction, but are not used in the calculation of the limit on 
a{pp -> if ) X Br{H — > 77). 

2. Correlated Systematic Uncertainties for CDF 

The dominant systematic uncertainties for the CDF analysis are shown in Table|lTl Each source induces a correlated 
uncertainty across all CDF channels' signal and background contributions which are sensitive to that source. 



3. Correlated Systematic Uncertainties for DO 
The dominant systematic uncertainties for the DO analysis are shown in Table IIIII 

TABLE IL Systematic uncertainties on the signal contributions for CDF's _H" — >■ 77 channels. Systematic uncertainties are listed 
by name; see the original references for a detailed explanation of their meaning and on how they are derived. Uncertainties are 
relative, in percent, and are symmetric unless otherwise indicated. 

CDF: // — >■ 77 channel relative uncertainties (%) 



Channel 



cc 


CP 


CC Conv 


CP Conv 




Systematic 


Uncertainties on Signal (%) 




6 


6 


6 


6 


14/7/5 


14/7/5 


14/7/5 


14/7/5 


2 


2 


2 


2 


3 


4 


2 


5 


3 


4 


2 


5 


0.2 


0.8 


0.1 


0.8 


- 


- 


0.1 


0.4 


0.2 


0.2 


0.2 


0.2 


- 


- 


7 


7 


0.4 


3.0 


0.2 


3.0 


1.0 


2.8 


1.0 


2.6 


3.0 


2.5 


1.5 


2.0 


0.4 


0.8 


1.5 


2.0 



Luminosity 

PDF 
ISR 
FSR 

Energy Scale 
Trigger Efficiency 
z Vertex 
Conversion ID 
Detector Material 
Photon/Electron ID 
Run Dependence 
Data/MC Fits 



Fit Function 



3.5 



Systematic Uncertainties on Background (%) 
1.1 7.5 



3.5 



TABLE III: Systematic uncertainties on the signal and background contributions for DO's _ff — > 77 channel. Systematic 
uncertainties for the Higgs signal shown in this table are obtained for mu = 125 GeV/c^. Systematic uncertainties are listed 
by name; see the original references for a detailed explanation of their meaning and on how they are derived. Uncertainties are 
relative, in percent, and are symmetric unless otherwise indicated. 

DO: H — >■ 77 channel relative uncertainties (%) 



Contribution 



Background 



Signal 



Luminosity 

Acceptance 

electron ID efficiency 

electron track-match inefficiency 

Photon ID efficiency 

Photon energy scale 

Cross Section 

Background subtraction 



2 
10 
3 

2 
4 

8 



3 
1 
10 



IV. COMBINED RESULTS 

Using the combination procedures outlined in Section III, we extract limits on SM Higgs boson production a x 
B{H — T' 77) in pp collisions at -^s = 1.96 TeV for 100 < rti}! < 150 GeV/c^. To facilitate comparisons with the 
SM and to accommodate analyses with different degrees of sensitivity, we present our results in terms of the ratio of 
obtained limits to the SM Higgs boson production cross section, as a function of Higgs boson mass, for test masses 
arranged in 5 GeV/c^ steps between 100 and 150 GeV/c^. 

The ratios of the 95% C.L. expected and observed limit to the SM cross section are shown in Figure [1] for the 
combined CDF and DO analyses. The observed and median expected ratios computed using the Bayesian and the 
CLg methods are listed for the tested Higgs boson masses in Table HVl In the following summary we quote only the 
limits obtained with the Bayesian method, which was chosen a priori. We obtain the observed (expected) values of 
7.2 (10.6) at TTiH = 100 GeV/c^, 10.5 (8.5) at m^ = 115 GeV/c^, and 13.4 (8.3) at mn = 120 GeV/c^. 

We also compute limits on the inclusive production cross section times the decay branching ratio cr(pp -> -ff + 
X) X Br{H —^ 77), in order to reduce model dependence from the SM predictions of the cross section and the decay 
branching ratio. Some residual model dependence remains, however, as the SM ratios of WH, ZH, VBF, and gg ^ H 
production are assumed in order to compute the signal acceptances for each channel. The signal acceptances vary 
by no more than a relative 10% between the different production mechanisms, however, and so the sensitivity to the 
predicted ratios is small. We perform the combination calculation for this result without the theoretical uncertainties 
on the total production cross section and the decay branching fraction Br{H — > 77). The resulting limits are shown 
in Figure [5] and listed in Table IVl 

TABLE IV: Ratios of median expected and observed 95% C.L. limit to the SM cross section prediction for the combined CDF 
and DO analyses as a function of the Higgs boson mass in GeV/c^, obtained with the Bayesian and CLs methods. 



Bayesian 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


Expected 
Observed 


10.6 
7.2 


9.2 
11.7 


8.7 
7.6 


8.5 
10.5 


8.3 
13.4 


8.2 
11.3 


8.9 
13.3 


9.6 
11.0 


10.8 
11.1 


12.7 
10.4 


16.4 
13.0 


CL. 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


Expected: 
Observed: 


10.5 
7.0 


9.3 

11.7 


8.9 
7.6 


8.5 
10.5 


8.2 
13.3 


8.2 
11.3 


8.7 
13.5 


9.4 
11.1 


10.7 
11.2 


12.9 
10.2 


16.7 
12.7 



TABLE V: The 95% CL limit on the product of the inclusive Higgs boson production cross section and the decay branching 
ratio to a pair of photons u{pp -^ H + X) x Br(H — > 77), as a function of the Higgs boson mass, assuming SM ratios of the 
four production mechanisms. 



Bayesian 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


150 


Expected (fb) 
Observed (fb) 


40.1 
27.2 


33.7 
42.6 


31.7 
26.8 


28.2 
34.9 


25.5 
42.0 


23.2 
31.9 


21.5 
32.6 


19.5 
22.5 


18.2 
18.6 


16.4 
13.6 


15.6 
12.4 



In summary, we combine CDF and DO results on SM Higgs boson searches with H — > 77, based on 7.0 fb^^ of data 
from CDF and 8.2 fb~^ of data from DO. We use the recommendation of the PDF4LHC working group for the central 
value of the parton distribution functions and uncertainties [2J1 . We use the highest-order calculations available for 
the gg — >■ iJ, WH , ZH , and VBF theoretical cross sections when comparing our limits to the SM predictions. We 
include consensus estimates of the theoretical uncertainties on these production cross sections and the decay branching 
fractions in the computations of our limits. 

The 95% C.L. upper limit on Higgs boson production is a factor of 10.5 times the SM cross section for a Higgs 
boson mass of mn =115 GeV/c^. Based on simulation, the corresponding median expected upper limit is 8.5 times 
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the SM cross section. Standard Model branching ratios, calculated as functions of the Higgs boson mass, are assumed. 
The results presented here extend significantly the sensitivity of the separate CDF and DO results. 
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FIG. 1: Observed and expected (median, for the background-only hypothesis) 95% C.L. upper limits on the ratios to the SM 
cross section, as functions of the Higgs boson mass for the combined CDF and DO analyses. The limits are expressed as a 
multiple of the SM prediction for test masses (every 5 GeV/c'^) for which both experiments have performed dedicated searches 
in different channels. The points are joined by straight lines for better readability. The bands indicate the 68% and 95% 
probability regions where the limits can fluctuate, in the absence of signal. The limits displayed in this figure are obtained with 
the Bayesian calculation. 
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FIG. 2: Limits on the product of the inclusive Higgs boson production cross section and the decay branching ratio to a pair 
of photons a-{pp -^ H -y- X) x Br{H — > 77), as a function of the Higgs boson mass, assuming SM ratios of the four production 
mechanisms. The bands indicate the 68% and 95% probability regions where the hmits can fluctuate, in the absence of signal. 
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